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Abstract 

ifj' production is studied in Pb-Pb collisions at 158 GeV/c per nucleon in- 
cident momentum. Absolute cross-sections are measured and production rates 
are investigated as a function of the centrality of the collision. The results are 
compared with those obtained for lighter colliding systems and also for the J 
meson produced under identical conditions. 
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1 Introduction 



The production of charmonium states in ultrarelativistic heavy ion collisions has been 
investigated since the mid 80s by the NA38 Collaboration [lj, followed in the 90s by 
NA50 [2]. The motivation is to search for the phase transition of nuclear matter 
from its normal state to a deconfmed quark-gluon plasma, as predicted to occur, un- 
der extreme energy densities or temperatures, by non-perturbative QCD. In such a 
deconfmed matter, charmonia states, in particular the J/0 vector- meson, have been 
predicted to be suppressed by Debye colour screening [3]. Indeed, the NA50 experi- 
ment has shown that J/0 production in Pb-Pb collisions is significantly suppressed 
with respect to the expectations derived from the rates measured in proton, oxygen 
and sulphur-induced reactions [U EJ E] • The 0' is a more loosely bound state than 
the J/0, which opens the possibility that it may dissociate at lower energy densities 
in connection with critical scenarios, as predicted by recent lattice calculations [3 [8], 
but also be sensitive to other mechanisms, as absorption in nuclear matter or disso- 
ciation by co- moving hadrons. After the first results on -0' production obtained by 
experiments NA38 [9] and NA50 [10], the framework for a general understanding of 
charmonium dissociation due to critical transitions has recently opened up with new 
experimental results on J/0 suppression [HI [12]. The measurement of ■0' produc- 
tion reported hereafter, performed with identical analysis procedures as for the 
should help and internally constrain the various theoretical approaches which try to 
give a coherent overall picture of charmonium production and suppression. 



2 Experimental setup, data selection and analysis 
method 

NA50 is a fixed target experiment at the CERN-SPS optimized for the detection of 
muon pairs produced in nucleus-nucleus collisions. The main component of the ap- 
paratus is a muon spectrometer made of an air-core toroidal magnet surrounded by 
two sets of multi-wire proportional chambers and trigger hodoscopes. An appropri- 
ate hadron absorber separates the spectrometer from the target itself, allowing the 
experiment to run with a high intensity incident beam and a moderate spectrometer 
illumination. Three independent detectors, located immediately downstream from 
the target region, measure the centrality of the reactions — a silicon strip multiplic- 
ity detector sampling the produced secondary charged particles, a Pb-scintillating 
fibers electromagnetic calorimeter measuring the neutral transverse energy Et and a 
very forward hadronic calorimeter (zero degree calorimeter) , embedded in the hadron 
absorber, which essentially measures the energy of the beam spectator nucleons in the 
collision, E Z dc- A beam hodoscope (BH), located 25m upstream from the target, 
allows to monitor and identify the incident beam and also provides an accurate time 
reference for trigger purposes. A detailed description of the NA50 apparatus can be 
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found in Ref. [T3] . 

This article is based only on the highest quality data samples collected by the 
experiment in years 1998 and 2000 when only one single target was used which, 
moreover, was placed in vacuum in year 2000. These features have efficiently solved 
problems affecting previous data collections such as secondary fragment reinteractions 
within the target assembly leading to a centrality smearing on one hand and Pb-air 
interactions simulating peripheral Pb-Pb collisions on the other. Furthermore, the 
selected data have benefited from an upgraded track reconstruction treatment based 
on improved software algorithms [H]. The data were collected with a typical 
beam intensity of 1-I.4xl0 7 ions/s, over 4.8 s bursts every 20 s, at 158 GeV per nu- 
cleon. The single Pb target was 4 mm thick in 2000 (3 mm in 1998), equivalent to 
10% (7%) of an interaction length. Muon pairs are selected in the rapidity window 
2.92 < yiab < 3.92 (0 < i/cm < 1) and with a Collins-Soper angle | cos9cs |< 0.5, 
resulting in an acceptance of around 14%. On-target interactions are selected requir- 
ing the proper correlation between hits in the two planes of the multiplicity detector. 
This target selection procedure is 88% efficient for E T > 5 GeV and 100% efficient 
for Et > 22 GeV. It allows to unambiguously identify the most peripheral on-target 
interactions and, therefore, has now been applied to the Pb-Pb data samples collected 
both in 1998 and in 2000. Moreover, muons produced off-target are further excluded 
by a cut on the transverse distance, in the middle target plane, between the muon 
track extrapolated from the spectrometer and the beam axis. Parasitic interactions 
occurring upstream from the target, mostly in the beam hodoscope, are rejected by 
a BH interaction detector and by anti-halo counters. Multiple piled-up interactions 
are discarded by a shape analysis of the signal from the electromagnetic calorimeter. 
Residual piled-up events are rejected by a 2cr cut in the Et — Ezdc correlation [2]. 

The opposite-sign dimuon invariant mass spectrum in the mass region of interest, 
as shown in Fig. [TJ results from five different contributions — the resonances J/ip 
and if)', as well as the continuum formed by Drell-Yan, by open charm semi-leptonic 
decays and by the unphysical combinatorial background originating mostly from un- 
correlated tt and K decays. In order to disentangle the different contributions and 
to provide the necessary ingredients for the determination of cross-sections, a Monte- 
Carlo simulation technique is used. It allows to determine both the acceptance of 
each physical process and the smearing of the corresponding detected muon pairs as 
they result from the actual experimental detector and selection criteria imposed on 
the data. Each physical process is separately generated and its corresponding muon 
pairs are propagated through the detector exactly the way real muons do, which ac- 
counts, in particular, for multiple scattering and energy loss in the absorbers and 
for the limited acceptance of the apparatus. The resulting smearing is particularly 
visible in the J/ip and ip' reconstructed shapes, entirely due to the experimental res- 
olution of the apparatus. The resonances are generated in rapidity with a Gaussian 
distribution of 0.6 units width. The transverse momentum used in the simulation 
follows the parametrization of a thermal distribution, dcr/dp T = p T M T K 1 (M T /T), 
where K\ is the modified Bessel function of the second kind and first order, in the 
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transverse mass, with T=236 MeV. These parameters were tuned in order to repro- 
duce the experimental spectra of the S-U collision system at 200 GeV/c [T5] and they 
account very well for Pb-Pb interactions at 158 GeV/c [16]. The Drell-Yan is evalu- 
ated with PYTHIA [17] , using the GRV LO 94 [18] set of parton density functions, 
with a minimum 4-momentum transfer and x domain as required by the NA50 phase 
space window. At the Born level, dimuons are produced with zero p?', to overcome 
this inadequate physical picture, a primordial kx inside the hadrons is considered, 
according to a gaussian with 0.8 GeV/c width. The DD mesons are also generated 
with PYTHIA, at leading order of QCD, using m c = 1.35 GeV/c 2 and a gaussian 
distributed with 1.0 GeV/c width [19] . The combinatorial background, free from 
any simulation, is completely determined both in shape and in amplitude from the 
measured like-sign pair distributions according to N^q = 2 V N ++ N . This method 
holds as long as the spectrometer acceptance does not depend on the muon charge, 
which is ensured by an appropriate selection of events. 

After the Monte-Carlo determination of the muon pair mass shapes of each phys- 
ical process, their corresponding amplitudes are obtained from a fit to the opposite 
sign dimuon invariant mass spectrum of the selected events. The background con- 
tribution, already established through the procedure explained above, is fixed. The 
DD contribution, although tiny in the mass range of relevance for the ip' amplitude 
determination and, moreover, practically uncorrelated with it, is roughly estimated 
from a separate fit in the dimuon mass range 1.7-2.2 GeV/c 2 . For this purpose, 
the corresponding Drell-Yan is deduced from the extrapolation of the muon pairs 
sampled above 4.2 GeV/c 2 , a mass region where Drell-Yan is the only contribution 
and is, therefore, unambiguously determined. The i/j', J/ip and the whole Drell Yan 
amplitudes are obtained from a fit in the mass range 2.9-7.0 GeV/c 2 , with both the 
combinatorial background and DD contributions fixed. The set of parameters which 
maximize the log-likelihood function is found by the MINUIT package [20] . which 
also evaluates the function in the neighbourhood of the minimum, thus providing the 
parameter statistical uncertainties. One should notice the good \ 2 °f the fit shown 
in Fig. [0 Details on the analysis method are given in [16] . 



3 Absolute cross-sections 

The study of ip' production in Pb-Pb collisions implies special care with the data 
treatment, in particular in what concerns the identification of systematic sources, 
since this resonance yields a weak signal due both to its small dimuon branching 
ratio {B^+^- = (7.3 ± 0.8) x 10~ 3 [2T] ) and to its large suppression with increasing 
centrality of the collision. Another challenge is the extraction of the ip' signal itself 
from the several overlapping dimuon sources. 

In order to measure absolute cross-sections, data collected in the year 2000 were 
selected under strict criteria with respect to the stability of the experimental condi- 
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Figure 1: The fitted opposite sign dimuon invariant mass spectrum, for mid-centrality 
Pb-Pb collisions. 

tions. After the data selection described in the previous section and after correcting 
for the target identification inefficiency, about 900 ip' are left for further analysis. 
This absolute number of events must still be corrected for various detector ineffi- 
ciencies, selection cut losses and acceptances in the NA50 phase space window, as 
detailed in Table [TJ The main systematical uncertainty contributing to the errors of 
the acceptance values is the generation model assumed. In particular, if a Fermi mo- 
tion model of the nucleons inside the nucleus is taken into account, the spectrometer 
acceptance of the resonances changes by ~1.5%. In the case of Drell-Yan, different 
parton density functions lead to a 4% difference in the acceptances. 

Absolute cross-sections for and J/i/j in the dimuon channel, as well as for Drell- 
Yan in the mass windows of 2.9 < M w < 4.5 and 4.2 < M w < 7.0 GeV/c 2 , are 
displayed in Table [2] within the kinematical ranges defined in section 2. The if}' ab- 
solute cross-section is affected by a systematic uncertainty of 7%, obtained as the 
quadratic sum of all uncertainties considered, namely the errors on the absolute nor- 
malizations, fit method, spectrometer acceptances and uncertainties associated with 
the Monte-Carlo inputs, like parton density functions and c-quark mass. The lat- 
ter, although directly related to Drell-Yan and open charm theoretical uncertainties, 
indirectly influence the ip' fitted normalization. The main contributions to the un- 
certainty in the ip' yield are the luminosity calculation, 4%, and the chosen PDFs to 
generate Drell-Yan, 5.8%, estimated as the numerical change in the result when using 
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Detector efficiencies (%) 


Dimuon trigger 


87 ± 1 


Muon track reconstruction 


94 ± 1 


Beam hodoscope 


83 ± 1 


Lifetime of DAQ 


96 ± 1 


Target identification (5< E T <22 GeV) 87 ± 2 


(Et >22 GeV) « 100 


Selection cut losses (%) 


Interaction pile-up 


15.8 ± 1 


Et — Ezdc correlation 


16.3 ± 1 


BH interaction 


2.2 ± 0.5 


Acceptances (%) 


if)' 


14.8 ± 0.3 


J/V 


12.5 ± 0.2 


Drell-Yan(2.9 < M w < 4.5) 


13.8 ± 0.2 


Drell-Yan(4.2 < M w < 7.0) 


17.8 ± 0.7 



Table 1: Detector efficiencies, selection cut losses and integrated acceptances. 
Cross-sections (fih) in < ycM < 1 an d | cos 9cs |< 0.5 
Pb + Pb — ► V' + x = °- 136 ± °- 013 -o°om 
Pb + Pb — ► J/ip + X = 23.32 ± 0.16 ± 0.84 
Pb + Pb — ► Drell-Yan(2.9 < M w < 4.5) + X = 1.37 ± 0.07 ± 0.05 
Pb + Pb — ► Drell-Yan(4.2 < M w < 7.0) + X = 0.172 ± 0.009 ± 0.009 

Table 2: Inclusive cross-sections for if)' and J /if), multiplied by their branching ratios 
into and for Drell-Yan, in the 4.2-7.0 GeV/c 2 mass range. 

CTEQ4L [22] instead of GRVL094. On the other hand, it has been verified that the 
more recent GRVL098 [23J induces a variation of less than 1.5% on the cross-section 
measurement. The systematic uncertainty due to the fit method is 2.8%. The J /if) 
cross-section obtained here, 23.32 ± 0.16 ± 0.84 fib, is 6.5% higher than the one 
derived from the Pb-Pb data collected in year 1995, 21.9 ± 0.2 ± 1.6 fib, reported in 
Ref. [13J. The difference results from the cumulative effects of the new J /if) experi- 
mental line shape [H] and the higher number of tracks reconstructed. Although the 
better description of the resonance tails would lead to a 3.5% smaller cross-section, 
the new offline software algorithms [5] deal much better with the high MWPC occu- 
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pancy levels induced by the Pb-Pb interactions and result in a net increase of the 
final cross-section value. 

Figure [2] shows the K exp factor, the ratio between the measured Drell-Yan cross- 
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Figure 2: The Drell-Yan K exp factor measured at 158-400 GeV/c beam momenta in 
p-A and A-B collisions. GRV L094 PDFs are used in the theoretical calculations 
and in the extraction of the measured values. 

section and the lowest order theoretical Drell-Yan cross-section. The constant be- 
haviour of this ratio from p-Be to Pb-Pb collisions shows that the Drell-Yan process 
is proportional to the number of nucleon-nucleon collisions in the NA38/NA50 phase 
space window. 

4 Cross-section ratios versus centrality 

The study of the ip' suppression as a function of centrality is based on 380 and 905 
ij}' events collected in 1998 and 2000, respectively, and analyzed as a function of their 
transverse energy used as the centrality estimator. Preliminary results of these stud- 
ies have already been presented in ref. [TUj . 

For each centrality class separately, the analysis is performed as described in sec- 
tion [21 In order to get the production pattern as a function of centrality, results in 
different classes have to be normalized with respect to each other. We have chosen 
to make use of the number of elementary nucleon-nucleon collisions for this normal- 
ization purpose and have replaced it by a proportional experimental direct estimate, 
namely the Drell-Yan cross-section in the corresponding centrality bin as determined 
from the same fit procedure. The method thus increases the robustness of the mea- 
sured centrality pattern, as corrections and efficiencies depending on centrality, if any, 
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cancel out in the ratio. 

Figure [3] displays the ratio B' a [if)') /a (DY) as a function of the transverse en- 




Figure 3: S^ctOOMDY) as a function of E T for the Pb-Pb 1998 and 2000 data samples 
(left). The cross-section ratios after combining the results obtained separately (right). 
Errors are the quadratic sum of statistical and systematic uncertainties. 

ergy, for the two data samples. We should note here that, although the target region 
was not in vaccum during the 1998 data taking period, the use of the multiplicity de- 
tector to tag in-target interactions has now allowed the study of peripheral collisions 
without ambiguities and, moreover, the results obtained from the year 1998 data 
sample are perfectly consistent with those obtained from the data collected in year 
2000, when the target was in vacuum [16]. Thus, the 1998 data have been rebinned to 
match the same 2000 centrality intervals in order to perform their weighted average. 
The combined result shows a suppression of the if)' production rate which becomes 
stronger with increasing centrality, up to a factor of 6 between peripheral and central 
collisions. Numerical values are reported in Table [3] Statistical and systematic errors 
have been quadratically added. PDF's used to generate Drell-Yan in the Monte Carlo 
simulations are the main contribution to the systematic errors, as explained in sec- 
tion [31 They are centrality dependent, but small when compared with the statistical 
errors. The Drell-Yan mass range chosen to normalize the if>' yield is 4.2-7.0 GeV/c 2 . 
One can scale up the results to the mass range used in the J/ip suppression studies, 
2.9-4.5 GeV/c 2 [2], applying the factor 7.96, which is the ratio between the Drell-Yan 
cross-sections in the two mass domains. 

The direct comparison between J /if) and if>' suppressions as a function of centrality 
is displayed on Fig. 0], which shows the ratio of cross-sections B'^aiif)') / B MM cr(J '/ 'if)) 
as a function of transverse energy. The corresponding numerical values are reported 
in Table [3j The statistical errors come mainly from the if)', since the J /if) statisti- 
cal uncertainty is less than 1%. Statistical and systematic errors have been added 
quadratically. From the steady decrease pattern of this ratio, we conclude that with 
respect to the J /if), the if)' is more and more suppressed with increasing centrality, 
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with a factor 2.5 between peripheral and central reactions. 

Table H] reports the number of ip', J ftp and Drell-Yan fitted events in each trans- 
verse energy interval. The numbers refer to both Pb-Pb 1998 and 2000 data taking 
periods. 
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Figure 4: B'^a^/B^a^/ip) as a function of E T for the Pb-Pb 1998 and 2000 data 
samples, separately for each year (left) and after combining both data sets (right). Errors 
are the quadratic sum of statistical and systematic uncertainties. 



Et range (GeV) 


(E T ) (GeV) 


BXV>')MDY4. 2 - 7 .o) 


WWVW) xl ° 3 


3-20 


13.9 


1.91 ± 0.35 


9.30 ± 1.24 to'ix 


20-35 


28.2 


1.48 ± 0.25 i£o4 


8.10 ± 1.32 + u f 9 


35-50 


43.0 


0.98 ± 0.18 


7.33 ± 1.35 t u f 7 


50-65 


57.8 


0.74 ± 0.17 t^i 


5.48 ± 1.36 t" f 3 


65-80 


72.7 


0.56 ± 0.15 


4.77 ± 1.40 ±^ 


80-95 


87.6 


0.59 ± 0.15 


5.59 ± 1.53 +^2 


95-150 


105.9 


0.31 ± 0.12 


3.56 ± 1.52 



Table 3: B'a(ijj') /UY 4.2-7.0 and B' a(ij)') / B^a^J /ip) values for each centrality bin, 
for the weighted average between the 2000 and 1998 results. The first error is statis- 
tical, the second systematic. 
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Fit ranee (OeV") 


ib' 


J /ib 


i-j J- 4.2 — 7. U 


3-20 


186 ± 25 


16942 ± 146 


112 ± 8 


20-35 


243 ± 31 


25229 ± 181 


187 ± 11 


35-50 


227 ± 35 


27276 ± 192 


264 ± 12 


50-65 


193 ± 36 


27681 ± 196 


288 ± 13 


65-80 


154 ± 36 


27315 ± 200 


310 ± 14 


80-95 


159 ± 37 


25111 ± 193 


311 ± 14 


95-150 


110 ± 40 


28570 ± 209 


401 ± 15 



Table 4: The numbers of ib', J/ifi and DY 4 . 2 _ 7 fitted events for each centrality range. 
Errors are only statistical. 

5 Comparison with ip' production in lighter colli- 
sion systems 

ip' suppression in Pb-Pb collisions can be compared with results obtained both from 
NA50 proton-induced and NA38 S-U reactions, the latter reanalysed by the NA50 
collaboration [24]. Data were collected with incident momentum beams of 400 and 
450 GeV/c for protons and of 200 AGeV/c for sulphur ions. In the framework of 
NA50, a coherent comparison can be done, since these data were collected with the 
same spectrometer, so that systematic effects related to trigger and reconstruction 
efficiencies are similar for all data sets. Besides, the transverse energy released in S-U 
collisions has been measured with a very similar electromagnetic calorimeter. Correc- 
tion factors have been applied to scale down the p-A and S-U center of mass energies 
and rapidity domains to the Pb-Pb conditions [25]. Taking into account the relative 
amount of protons and neutrons in each colliding nucleus, isospin corrections are 
needed when the Drell-Yan cross-section is compared among different systems. For 
this reason the measured Drell-Yan cross-sections have been normalized to proton- 
proton collisions. One should notice that, since ip' production in the NA38/NA50 xf 
domain is dominated by gluon fusion, rather than by quark-antiquark annihilation, 
the ib' isospin corrections are negligible. 

In order to compare the ip' behaviour between different colliding systems we con- 
sider hereafter three variables directly related to the impact parameter of the inter- 
actions, namely L, the average path crossed by the cc pair inside the nucleus, N part , 
the number of nucleons participating in the interaction, and e, the reached energy 
density. They have been calculated through the Glauber model formalism [26] of 
nucleus-nucleus collisions. Such a model considers an A-B interaction as a superposi- 
tion of independent interactions between the beam and target nucleons. The model 
also assumes that the nucleon-nucleon cross section remains unchanged as the two nu- 
clei cross each other. One should bear in mind that this is an approximation because 
a nucleon, after a collision, may become excited and subsequently interact with other 
nucleon with a different cross section. The Pb nuclear density is parameterised by 
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a two-parameters Fermi function [27], also known as the Woods-Saxon distribution, 
p(r) = Po/(l + exp((r — r )/c)), where p is the average nuclear density, taken as 
0.17 fm -3 , ro is the half-density radius of 6.624 fm and c is a diffuseness parameter 
of 0.549 fm for protons. Taken into account the "neutron halo" effect [28j [6] a dif- 
ferent diffuseness parameter (0.667 fm) is used to describe the neutrons distribution 
inside the Pb nucleus. The same model has been applied to the U nucleus, with a 
half-density radius of 6.8054 fm and diffuseness parameters of 0.605 fm for protons 
and 0.786 fm for neutrons. The actual shape of the U nucleus is not taken into ac- 
count. Figure [5] shows the ratio B'^aiip') / 'cx(DY) as a function of L. The measured 
suppression patterns suggest the following three features: a) a fair agreement with 
exponential behaviours; b) two different regimes, one for proton and a different one 
for ion-induced reactions; c) a similar centrality dependence for S-U and Pb-Pb inter- 
ations. Using an exponential parametrization to describe if)' absorption as a function 
of L [29], according to exp(— (pL)cr a b s ), the fit of the data gives an absorption cross- 
section of 7.3 ±1.6 mb in p-A collisions, while a much higher value, 19.2 ± 2.4 mb, 
is obtained for ion- ion collisions (S-U and Pb-Pb fitted simultaneously). The left 




Figure 5: B'a(ip') / 'cr(DY) as a function of L. Statistical and systematic uncertainties are 
added quadratically. 



panel of Fig. O shows the same ratio as a function of the number of participants. By 
definition, a nucleon is designated as participant if it has at least one inelastic inter- 
action with one or more surrounding nucleons. For this work the calculation of N part 
has followed eq. 9 of Ref. [30]. It is worth mentioning that, at SPS energies, several 
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experiments have observed that the number of participating nucleons in a nucleus- 
nucleus collision scales with the transverse energy [311 The ip' suppression as a 
function of the energy density is shown in the right panel of Fig. [6j but only for A-B 
collisions. In the framework of the Bjorken model [33], the energy density reached 
in a nucleus-nucleus collision is proportional to the number of participating nucleons 
per unit transverse area (see Refs. [3U [TH] for examples of such an evaluation). An 
initial formation time, r = 1 fm/c, has been assumed. It is important to keep in 
mind that the three centrality variables, L, N part and e, are strongly correlated, 
since they are all calculated by the same Glauber model of nucleus-nucleus collisions. 
Their average and r.m.s. values for each centrality bin are quoted in Tables [5] and [61 
for Pb-Pb and S-U colliding systems, respectively. These values take into account the 
smearing due to the experimental resolution of the electromagnetic calorimeters used 
by NA50 and NA38. Concerning the energy density values, we are quite confident 
that the comparison between S-U and Pb-Pb is robust, while the absolute values 
could suffer from a 20-25% systematic uncertainty. 

In principle, the use of several centrality estimators in the study of the ip' suppres- 
sion pattern in two significantly different colliding systems, such as S-U and Pb-Pb, 
should indicate the centrality variable as a function of which the two measured sup- 
pression patterns would best overlap. In practice, however, the statistical errors of 
the ijj' suppression patterns presently available do not allow us to make very clear 
statements in this respect, apart from saying that the use of the N part variable seems 
to result in a somewhat worse overlap between the two data sets. The Pb-Pb and S-U 
slopes, from a linear fit to the ip'/UY ratio vs. N part or e in the overlapping centrality 
range, differ by 1.6 (resp. 1.0) standard deviations. These observations could also be 
slightly affected by the approximate nature of our calculations in what concerns the 
shape of the U nucleus, as mentioned before. 



Et range 


{E T ) 


(L) 


r.m.s 


\Npart) 


r.m.s. 


<*> 


r.m.s. 


(GeV) 


(GeV) 


(fm) 


(fm) 






(GeV/fm 3 ) 


(GeV/fm 3 ) 


3-20 


13.9 


4.90 


0.84 


44.6 


17.7 


1.24 


0.37 


20-35 


28.2 


6.65 


0.44 


96.7 


18.0 


2.04 


0.20 


35-50 


43.0 


7.65 


0.31 


147.1 


19.3 


2.53 


0.14 


50-65 


57.8 


8.34 


0.24 


197.7 


20.6 


2.89 


0.11 


65-80 


72.7 


8.83 


0.18 


248.6 


22.0 


3.19 


0.09 


80-95 


87.6 


9.17 


0.14 


299.2 


23.0 


3.52 


0.07 


95-150 


105.9 


9.43 


0.10 


353.3 


22.4 


3.76 


0.06 



Table 5: Mean free path crossed by the cc pair inside the nucleus, number of parti- 
cipating nucleons and energy densities for each Pb-Pb Et range. 
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~p]rn VP\ n P*P 




(L) 


r.ixi. s 


(N A 


r.m.s. 


(e) 


r.m.s. 


(GeV) 


(GeV) 


(fm) 

V 1111 / 


(fm) 






(GeV/fm 3 ) 

\ * — ' Vj v / nil / 


(GeV/fm 3 ) 

\ v_ * \_, v / 1111 / 


13-28 


22.1 


4.37 


0.54 


31.3 


9.1 


1.04 


0.22 


28-40 


34.3 


5.35 


0.43 


50.0 


9.7 


1.46 


0.18 


40-52 


46.3 


6.04 


0.41 


66.7 


10.8 


1.76 


0.16 


52-64 


58.3 


6.65 


0.41 


83.5 


11.7 


2.01 


0.15 


64-76 


70.2 


7.17 


0.38 


98.8 


11.0 


2.22 


0.13 


76-88 


81.7 


7.53 


0.29 


109.0 


8.2 


2.38 


0.09 



Table 6: Same as previous table, for S-U collisions. 
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Figure 6: B'a(ip') /cr(DY) as functions of the number of participating nucleons (left) and of 
the energy density (right). Statistical and systematic uncertainties are added quadratically. 



6 Comparison between J/ip and ip f suppressions 

The ratio B'a(iff) / 'S MAt cr(J /ip) as plotted in Fig. [7] shows the relative suppression of 
the two charmonium states as a function of the product of the projectile and target 
mass numbers. The ratio is parametrized, in the range of the p-A data points, by the 
power law A Aa , where a accounts for all nuclear effects suffered by the resonances. 
The measured negative value of this parameter, Aa = —0.048 ± 0.008, shows that 
the ij}' is more absorbed than the J/ip, already in p-A collisions [5]. Furthermore, a 
strong suppression of the ip' with respect to the J/ip is measured for the ion induced 
reactions, S-U and Pb-Pb. 

Figure [8] shows, for various interacting systems, the ratio between the mea- 
sured charmonium yields, normalized to Drell-Yan dimuons, and the corresponding 
expected "normal nuclear absorptions". The latter are computed from the p-A re- 
sults. Using full Glauber calculations, the absorption cross-sections in nuclear matter 
are cr^ s b (J/^) = 4.2 ± 0.5 mb and cr^O') = 7.7 ± 0.9 mb, respectively HQ [6]. One 
should notice the good compatibility between this cr^i'ip') value and the one obtained 
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Figure 7: B'a(ip')/ 'B^trfJ /ip) as a function of A pro j ec tii e xBtarget- Statistical and system- 
atic uncertainties are added quadratically. 



with a simple exponential in L, 7.3 ± 1.6 mb, as reported in section [5j The figure 
shows that the departure from ordinary nuclear absorption occurs at a lower value of 
L for ip' than for J/ip. Moreover, while the ip'/UY behaviour changes between p-A 
and the whole centrality range of ion-ion interactions, the J/^/DY trend in Pb-Pb 
collisions begins to be compatible with both p-A and S-U results, but deviates more 
and more as the collision centrality increases. More details on the J/ip anomalous 
suppression in Pb-Pb interactions can be found in Ref. [2J. 



Conclusions 

We have measured the ip' absolute cross-section in Pb-Pb collisions at y/s = 17.3 GeV, 
in the NA50 kinematic domain, B'^a{^') = 0.136 ± 0.013(stat)+0.010 - 0.006(syst) 
fib. The study of the ip' production as a function of centrality shows a steady decrease 
from peripheral to central collisions. It leads to a ip' suppression factor of 6 relative to 
Drell-Yan production and of 2.5 relative to J/ip. The comparison with lighter systems 
shows that the ip' is much more suppressed in nucleus-nucleus than in proton-nucleus 
reactions and that the suppression pattern in S-U and Pb-Pb approximately overlap 
when studied as a function of L or energy density, while the overlap as a function of 
the number of participating nucleons seems to be somewhat worse. The comparison 
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Figure 8: The ratio "measured over expected" for the relative yields B' flfl a('ip')/a(DY) and 
B^^<7(J/ip)/a(DY), as a function of L. 



between the J/ip and ip' suppression patterns shows that the ip' anomalous suppres- 
sion sets in for lower values of L. 
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